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Results on a search for pair production of second generation scalar leptoquark in pp collisions 
at y^=1.96 TeV are reported. The data analyzed were collected by the CDF detector during the 
2002-2003 Tevatron Rnn II and correspond to an integrated luminosity of 198 pb“^. Leptoqnarks 
(LQ) are songht throngh their decay into (charged) leptons and quarks, with hnal state signatures 
represented by two muons and jets and one muon, large transverse missing energy and jets. We 
observe no evidence for LQ production and derive 95% C.L. npper limits on the LQ prodnetion 
cross sections as well as lower limits on their mass as a function of f), where jd is the branching 
fraction for LQ —> p,q. 

PACS numbers: 14.80.-j, 13.85.Rm 


The symmetry between leptons and quarks present in 
the standard model (SM) of particle physics, has lead 
to the theoretical speculation that a more foundamental 
force could operate at energy scales larger than the elec- 
troweak symmetry breaking scale, allowing for quark to 


lepton transitions, mediated by new gauge bosons. The¬ 
ories like grand uniheation or R-parity violating super- 
symmetric models introduce the idea of quark to lepton 
transitions^. Whenever quarks and leptons are allowed 
to couple directly to each other, new bosons carrying 


4 


both lepton and baryon quantum numbers can also ex¬ 
ist. They are called leptoquark (LQ)0, they can have 
spin 0 (scalar LQ) or 1 (vector LQ),are color-triplet parti¬ 
cles, and can either be produced singly or in pairs. Most 
of the other characteristics, such as weak ispospin, elec¬ 
tric charge and the coupling A to lepton and quark are 
model dependent 0. Their masses are not predicted. To 
accomodate experimental constraints on flavor changing 
neutral currents, LQ are assumed to couple to fermions 
of the same generation^ 0. While first generation LQ 
have been extensively searched for at e — p, e'^e~ and 
pp colliders, via single and pair production, and strong 
limits on their production cross section and mass have 
been set0, 0, second and third generation LQ are de¬ 
tectable via pair production at hadronic colliders and can 
be singly produced at e+e” machines. Indeed upper lim¬ 
its on second generation LQ production cross section and 
lower limits on their masses exist from LEP0 and the 
Tevatron Run 100, but they are usually weaker than 
the ones for first generation LQ. In this paper we present 
new results on a search for second generation LQ pair 
production in pp collisions at •\/s=l-96 TeV. Assuming 
two different values for the LQ branching fraction (3 to 
muon and quark, we consider the following final state sig¬ 
natures: both LQ decaying into muons of opposite charge 
and quarks (/3 = 1), and one of the LQ decaying into a 
muon and a quark and the other into a neutrino and a 
quark (/3 = 0.5). Since we do not observe evidence for 
LQ production, we set an upper limit on the production 
cross section times branching fraction (Br) where Br = 
(33 or Br = 2/3(1 — (3). These results are then combined 
with the one from a search for scalar LQ pairs decaying 
into WQ Q, re sulting in jets and missing transverse energy 
topology [lOj. In this way we can express our limits as a 
continous function of the parameter (3. 

CDF is a general-purpose detector built to study the 
physics of pp collisions at the Tevatron accelerator at Fer- 
milab and it is described in detail elsewhere . We use 
a cylindrical coordinate system around the beampipe in 
which 6 is the polar angle, (j) is the azimuthal angle and 
77 = —ln(tan(|)). The transverse energy is defined as 
Et = E sin 9 and the transverse momentum is defined 
as Pt = Psin0, where E is the energy measured by the 
calorimeter and P the momentum measured by the track¬ 
ing system. The vector I(t is equal to —EE^ni where 
is a unit vector that points from the interaction vertex 
to the ith calorimeter tower in the transverse plane. Its 
magnitude, 1(!t, is called transverse missing energy. 
is corrected following the correction of jet energies, and 
if muons are identified in the event, I^t is corrected for 
the muon momenta. 

The data used in the analysis were collected during 
the 2002-2003 Tevatron Run 11. The integrated luminos¬ 
ity for this data sample is 198 ± 12 pb“^. Events are 
selected online by requiring track segments (“stubs”) re¬ 
constructed in the muon chambers and matched to indi¬ 


vidual tracks reconstructed in the central tracker (high 
Pt muon trigger). The efficiency of the trigger combi¬ 
nations used in the uu ii a nd pvjj analyses, measured 
using Z ^ data[^ |^, is ^90%, varying from 

about 87% to 95% depending on the type of muon cham¬ 
ber used to detect the candidate muon. Muons are se¬ 
lected as “tight” or “loose” (this second category being 
used in the pp analysis only). A tight muon requires a 
reconstructed track segment in the muon chambers with 
positions well matched to the extrapolation of a single 
track, while a loose muon is the one selected by requiring 
only one isolated track. In both cases the energy deposi¬ 
tion in the calorimeters must be consistent with that of 
a minimum-ionizing particle. We apply a cut on the 
of the track fit to eliminate kaons and pions which have 
decayed in flight. The identification efficiency for muons 
has also been measured using data ^,0 and is approx¬ 
imately 90%, going from 89% to 95% for different types 
of detector and selection used to identify the candidate 
muon. The coordinate of the lepton (also assumed to be 
the event coordinate) along the beamline must fall within 
60 cm of the center of the detector ( Zyertex cut) to ensure 
a good energy measurement in the calorimeter. This cut 
has an efficiency of (95 ± O.I(stat) ± 0.5(sj/s))%, and is 
determined from studies of minimum bias events. The 
efficiencies of the identification cuts, the trigger selection 
and the vertex cut, measured using data are taken into 
account by using scale factors between data and Monte 
Carlo events. Jets are reconstructed using a cone of fixed 
radius R = (A? 7 ^ J- A(/)^) = 0.7 and for these analyses 

are required to be in the |? 7 | < 2.0 range. Jets are cali¬ 
brated as a function of 77 and Et and their energy is cor¬ 
rected to the parton level [10. Neutrinos produce missing 
transverse energy, IpT-, which is measured by balancing 
the calorimeter energy in the transverse plane. The muon 
sample is heavily contaminated by events produced by 
cosmic rays interactions with the detector. Since these 
events do not originate from a common interaction ver¬ 
tex, the timing capability of the Central Outer Tracker 
(COT) is used to reject events with two muon tracks, one 
of which travels toward the beam pipe. We also require 
that the muon track passes close to the beam line, within 
distances less than 0.02 cm ( 0.2 cm) for tracks with (with¬ 
out) silicon hits. In the analyses we are describing, the 
signal selection criteria are set according to the kinematic 
distribution (e.g. pt of the muons and Et of the jets) 
of decay products determined from Monte Carlo studies, 
optimized to eliminate background with a minimal loss 
of signal events [l^. 

In the dimuon J- jets topology, from the inclusive muon 
triggers dataset we select events with two reconstructed 
isolated muons with Pt > 25 GeV/c. The first muon 
is required to be tight, i.e. to have a stub associated 
to a track, while the second one can be without a stub 
(“stubless”). Events are further selected if there are at 
least two jets with Et > 30 and 15 GeV respectively. In 
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the search in the muon, neutrino and two jets topology, 
we select events with one reconstructed tight muon with 
Pt > 25 GeV/c. We veto events with a second loose or 
tight muon to be orthogonal to the previous selection. We 
then accept events where there is large missing transverse 
energy, Ij^T > 60 GeV and at least two jets with Et >30 
GeV. 

The above datasets are composed predominantly of 
events coming from QCD production of Z/W bosons in 
association with jets and tt production (where one or 
both the W's from top decay into muon and neutrino). 
To reduce these backgrounds we apply several cuts which 
depend on the final state topology. 

1 ) /ip, analysis: 

i) veto of events whose reconstructed dilepton mass 
falls in the window 76 < < 110 GeV/c^ to remove 

the Z + jets contribution and < 15 GeV/c^ to avoid 
contamination from J/di and T production; ii) Exiji) + 
Et{ 3 2 ) > 85 GeV and Eriei) + ET{e 2 ) > 85 GeV, 
hi) \/ {{Eriji) + ET{j 2 )Y + {Erie-i) + ET{e 2 )Y) > 200 
GeV. The effect of the last two cuts is shown in Fig. 1, 
where SM background is compared to a LQ signal for 
Mlq = 220 GeV/c 2 . 
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FIG. 1: Graphical representation of the last two topological 
cuts applied in the ppjj analysis as observed on MG events. 
The discrimination between SM background and LQ signal is 
evident. 


2 ) fj.iy analysis: 

— jet) > 5° to veto events where the trans¬ 
verse missing energy is mis-measured due to a mis- 
measure of the jet energy, and A(f){IfT — p) < 175° to 
ensure that that the missing energy does not come from 
mis-measurement of the muon momentum. ii)i?T(ji) + 
ET{j 2 ) > 80 GeV, iii)MT(pj^) > 120 GeV/c^ to reduce 
the W -b 2 jets background; iv) a mass dependent cut con¬ 
sisting in selecting events falling in mass windows defined 
around several LQ masses. We require that the recon¬ 
structed mass combinations of the jets, muon and I/^t be 
consistent with those reconstructed from the LQ Monte 
Garlo. The ambiguity of the jet assignements allows for 
two different sets of reconstructed LQ pairs in each event, 


from which we build lineshapes of mass distributions by 
matching the reconstructed objects to the generator level 
objects. The lineshapes are then used with the following 
selection: |M(p, ji)-MLQ| < 2cri or j 2 )- Mlq\ < 

2 a 2 , and Ji) > T™™ or J 2 ) > 1 ^“, 

where Mlq is the mean of the reconstructed LQ distribu¬ 
tion and cti ^2 are the width parameterizations. M(p, ^ 1 ^ 2 ) 
are the reconstructed muon-jet invariant masses while 
31 , 2 ) are the reconstructed neutrino-jet trans¬ 
verse invariant masses. For the transverse mass distri¬ 
butions, we have chosen a mass-dependent lower cut de¬ 
noted by TYjr- = 20 -b {Mlq - 120) GeV/c^ and 
r 2 ™" = 20 -b {Mlq - 120)/2 GeV/c^. In Fig. 2 we plot 
the mass distributions of the selected events (before the 
mass limit cut) compared to the signal distribution for 
rriLQ = 180 GeV/c^. 





FIG. 2: Final mass distributions (see text) of the surviv¬ 
ing events before the mass limit cut compared to the signal 
distribution for Mlq = 180 GeV/c^. 

We study the properties of the physics backgrounds by 
generatin g ev ents corresponding to Z/W + 2 jets with 
ALPGEN|16| -b HERWIG fljl f to perform parton show¬ 
ering) and tt with PYTHIA [ij . A complete simulation 
of the GDF II detector based on GEANt[i^ and full 
event reconstruction is then performed. To normalize the 
number of simulated events to data we use the theoretical 
cross sections for tt from and for 7 /Z ^ pp -b 2 jets 
from |2l| . The background arising from multijet events, 
where a jet is mismeasured as a muon or where the muon 
comes from pion decay (QCD/fake), is evaluated using 
data. In the ppj j analysis we examine the data for same- 
sign events (events with two muons of the same charge) 
remaining after each kinematical cut. We estimate the 
background contribution to be twice the number of same 
sign events, in the assumption that there is no evidence of 
LQ signal in these type of events (the LQ pair have oppo¬ 
site charge, giving rise to two opposite charge muons). In 
the pi^jj analysis the contribution from the QCD/fakes 
background is estimated by examining the phase space of 
the Ij^T vs the muon fractional isolation for data events 
in which the muon isolation requirement is not enforced. 
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Here the muon fractional isolation is defined as the ra¬ 
tio between the calorimetric energy not associated with 
the lepton in a cone of Ai? = 0.4 around the lepton and 
the energy of the lepton. The following assumptions are 
made: since jets are produced in association with other 
particles, the isolation fraction of a jet will generally be 
larger than the one corresponding to a muon; there is no 
correlation between the isolation of the muon and IpT, 
and in the region where the Ifr is small and the isola¬ 
tion of the muon is large the LQ contribution is expected 
to be negligible (background-dominated region). With 
these assumptions, from the ratio of the number of events 


in the background-dominated regions we can extrapolate 
the contribution in the signal region. Other backgrounds 
from bb, Z ^ rf, WW are negligible due to the muon 
isolation and large muon and jet transverse energy re¬ 
quirements. In the /i/i channel the expected number of 
Z 2 jets events is 1.7 ±0.1. The expected number of tt 
events is 0.22±0.03 events. We estimate 1±1 fake events 
The overall background estimate is: 3 ± 1 events. In the 
fj,iy channel, the number of events in each mass region, 
compared with the background expectations is reported 
in Table I. 


TABLE I: Number of events surviving all cuts in the muon, missing energy and jets topology, compared with background 
expectations, as a function of the LQ mass (in GeV/c^). Errors on the expectations are both statistical and systematic. 


Mass 

140 

160 

180 

200 

220 

Wjj 

0.9 ± 0.1 

1.4 ± 0.1 

1.4 ± 0.1 

1.6 ± 0.1 

1.6 ± 0.1 

top 

1.7 ± 0.2 

1.8 ± 0.2 

1.4 ± 0.2 

1.0 ± 0.2 

0.8 ± 0.2 

Zjj 

0.20 ± 0.01 

0.20 ± 0.01 

0.20 ± 0.01 

0.20 ± 0.01 

0.20 ± 0.01 

multijets 

0.3 ± 0.3 

0.3 ± 0.3 

0.3 ± 0.3 

0.3 ± 0.3 

0.3 ± 0.3 

Total 

3.1 ± 0.3 

3.7 ± 0.4 

3.2 ± 0.3 

3.1 ± 0.3 

2.9 ± 0.3 

Data 

3 

4 

2 

0 

0 


We check the prediction of our background sources with 
data in control regions where the background contribu¬ 
tion is maximized. For the analysis the region is 
defined by requiring two muons with Pt > 25 GeV/c, 
75 < < 105 GeV/c^ and 2 jets with Et > 30,15 

GeV. We observe 110 events and expect 88 ± 10. For the 
analysis we ask for one muon with Pt > 25 GeV/c, 


> 35 GeV and 2 jets with Et > 30 GeV and observe 
203 events to be compared with a prediction of 221 ± 15 
from SM sources. 

The efficiency to detect our signal is obtained from MG 
simulated LQ (PYTHIA) events to account for kinemat- 
ical and geometrical acceptance. The total efficiencies for 
a LQ signal are reported in Table 11. 


TABLE II: Efficiencies after all cuts, errors (statistical and systematic) and 95% C.L. upper limits on the production cross 
section x branching fraction Br, as a function of M^q, for the two channels. 


MLQ(GeV/c2) 

mj 





€ 

(TxBr(pb) 

e 

CTxBr(pb) 

100 

0.020 ± 0.003 

1.35 

0.0050 ± 0.0005 

- 

120 

0.05 ± 0.005 

0.52 

0.070 ± 0.005 

0.86 

160 

0.13 ± 0.01 

0.18 

0.070 ± 0.005 

0.73 

200 

0.19 ± 0.02 

0.13 

0.110 ± 0.005 

0.41 

220 

0.21 ± 0.02 

0.11 

0.13 ± 0.01 

0.24 

240 

0.24 ± 0.02 

0.10 

0.13 ± 0.01 

0.24 

260 

0.26 ± 0.02 

0.09 

0.14 ± 0.01 

0.21 


The following systematic uncertainties are considered when calculating signal acceptance and background pre- 
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dictions: i) luminosity: 6% ii) choice of parton distribu¬ 
tion functions: 2.1% iii) statistical error of MC < 1.2% 
iv) jet energy calibration scale < 1% v) muon reconstruc¬ 
tion : 0.8% vi) Zyertex cut : 0.5%. vii) initial and final 
state radiation 1.8%. After all selection cuts, 2 events 
remain in the /i/i channel, while the number of events 
remaining in the channel is reported in Table I. 

In the analyses described above the number of events 
passing the selection cuts is consistent with the expected 
number of background events. The conclusion of the two 
searches is that there is no LQ signal: hence we derive 


an upper limit on the LQ production cross section at 
95% confidence level. We use a Bayesian approach 
with a flat probability distribution for the signal cross 
section and Gaussian distributions for acceptance and 
background uncertainties. The cross section limits are 
tabulated in Table II and the mass limits are tabulated 
in Table III. To compare our experimental results with 
the theoretical expectation, we use the next-to-leading 
order (^LO) cross-section for scalar LQ pair production 
from [i^ with CTEQ6 parton distribution functionsjl^. 



FIG. 3: 95% C.L. limit on the experimental cross section times branching ratio as a function of the LQ mass for the i-ivjj and 
/r/rj’jchannel. The NLO theoretical cross section is plotted for different values of the renormalization scale. Mass limits of 170 
GeV/c^ and 224 GeV/c^ respectively are obtained. 


The theoretical uncertainties correspond to the variations 
from Mlq/2 to 2Mlq of the renormalization scale /i used 
in the NLO QCD calculation. To set a limit on the LQ 
mass we compare our 95% CL upper experimental limit 
to the theoretical cross section for /r = 2 Mlq, which is 
conservative as it corresponds to the lower value of the 
theoretical cross section. We find lower limits on M(LQ) 
at 224 GeV/c2 (/3 = 1) and 170 GeV/c^ (/3 = 0.5). They 
are reported in Figure 3. To obtain the best limit how¬ 
ever, we combine the results from the two decay channels 
just described with the result of a search for LQ in the 
case where the LQ pair decays to a neutrino and quark 
with branching ratio Br{LQ vq) =1.0^3- The in¬ 
dividual channel analyses are in fact optimized for fixed 
values of (3 (1,0.5,0) while in the combined analysis, due 
to the contributions of the different decay channels, the 
signal acceptance can be naturally expressed as a func¬ 
tion of /3. As for the treatment of uncertainties, the 
searches in the and pLvjj channels use common cri¬ 
teria and sometime apply the same kind of requirements 
so the uncertainties in the acceptances are considered 
correlated. When calculating the limit combination in¬ 
cluding the vujj channel the uncertainties are considered 


uncorrelated. For each /3 value a 95% C.L. upper limit 
on the expected number of events is returned for each 
mass, and by comparing this to the theoretical expecta¬ 
tion, lower limits on the LQ mass are set. The combined 
limit as a function of (3 is shown in Figure 4, together 
with the individual channel limits. The combined mass 
limits are also tabulated in Table III. 

The final result presented here is better than the results 
obtained with Tevatron Run I dataj^ . This is mostly due 
to the small increase in the cross section as a function of 
the center of mass energy (from 1.8 to 1.96 TeV), and 
an increase in the muon acceptance. A comment is in 
order when comparing this result with the ones recently 
published^ for first generation LQ and third generation 
lqIH. While the signatures of LQ production is very 
similar (high Pt leptons, large transverse missing energy 
and energetic jets) one has to consider the constraint 
on the LQ particle to only couple to same generation 
fermions. This implies different types of selection and 
exclude the possibility of combining inter-generation re¬ 
sults. Also, in the case of electrons and muons, as can be 
seen from the current results, the similarity in their ac¬ 
ceptance results in similar cross section and mass limits. 
























TABLE III: 95% C.L. lower limits on the second generation 
scalar LQ mass (in GeV/c^), as a function of /3.The limit from 
cdf[ 3 Run I (~ 120p6 is also given. 


p 

MM Jj 


vvjj 

Combined 

CDF Run I 

0.01 

- 

- 

114 

125 

- 

0.05 

- 

- 

110 

133 

- 

0.1 

- 

137 

- 

143 

- 

0.2 

- 

155 

- 

157 

- 

0.3 

100 

162 

- 

176 

- 

0.4 

152 

168 

- 

200 

- 

0.5 

171 

170 

- 

208 

- 

0.6 

184 

168 

- 

213 

- 

0.7 

196 

162 

- 

217 

- 

0.8 

206 

155 

- 

221 

- 

0.9 

215 

137 

- 

224 

- 

1.0 

224 

- 

- 

226 

202 


Search For Second Generation Scalar Leptoquarks 



Leptoquark Mass (GeV/c^) 


FIG. 4: Leptoquark mass exclusion regions at 95% C.L. as 
function of Br(LQ^ fj, q). 

while the result is quite different in the case of third gen¬ 
eration LQ, due to the much smaller acceptance for t lep¬ 
tons. The results presented in this paper, as well as the 
ones recently published in 0 , are obtained with a statis¬ 
tical sample corresponding to about twice the luminosity 
collected in Tevatron Run I. A future increase in the Run 
II luminosity by an order of magnitude is estimated to 
extend the LQ mass range to ^ 300 GeV/c^ for the /3 
= 1 case M- Substantially higher LQ masses will be 
explored at the future Large Hadron Collider (LHC)!!!. 

In conclusion, we have performed a search for pair pro¬ 
duction of second generation scalar LQ in the dimuons 
-I- jets and muon, missing energy -|- jets topologies, us¬ 
ing 198 pb“^ of proton-antiproton collision data recorded 
by the CDF experiment during Run II of the Tevatron. 


We combined these findings with the ones from a search 
in the I/^t + jets topology [ifij. No evidence for LQ is 
observed. Assuming that the LQ decays to muon and 
quark with variable branching ratio /3 we exclude LQ 
with masses below 226 GeV/c^ for /3 = 1, 208 GeV/c^ 
for P = 0.5 and 143 GeV/c^ for /3 = 0.1 at 95% C.L. 
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